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Mass of nucleus in units of the neutron mass. 


Act 


Activation ** rate of disintegration. 


o< 


Efficiency of counter in a particular geometric 
arrangement. 


B 


Arbitrary constant. 


(3 


Negative beta particles. 


E 


Energy of neutron. 


El 


Energy of neutron at birth. 


ev 


Electron volts. 


$(e,/) 


Neutron flux distribution = number of neutrons per 
3 

cur at t cm from the source with energy between E 
and E + dE. 


A 


Decay constant. 




Neutron density. 


N 


3 

Density of nuclei = number per cm . 


'n. 


Activation due to resonance neutrons divided by 


q(E,r) 


activation due to neutrons of all energies. 

3 

Slowing down density = numbers of neutrons per cnr 
at distance V cm from the source slowing down past 
energy E per sec. 


Qo 


Source strength = number of neutrons released 
per sec. 


r 


Distance from the source. 


(HL 


Absorption cross section. 


*7 


Scattering cross section. 
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I Introduction 



In a nuclear power reactor neutrons emitted from the fission process of uranium 
235 must be slowed down inside of the reactor by a moderator if they are to propagate 
the fission reaction. Because of the importance in reactor design of knowing the 
spatial distribution of these slowed neutrons many studies have been made of the 
slowing down process. Probably the earliest successful method of treating the slow- 
ing down theory can be attributed to Amaldi and Fermi. In their model* collisions 
between neutrons and the nuclei of the moderator are treated as free particle col- 
lisions o 

Bethe^^ has calculated that the effeot of chemical binding of the moderator 
atoms can be neglected in the treatment of the collisions between neutrons and modera- 
tor atoms. The Fermi slowing down theory should therefore be valid for the slowing 
down of neutrons whose energy is less than the binding energy of a moderator atom in 
its lattice. The purpose of this experiment is to test the validity of the Fermi 
theory when applied to neutrons whose energy is less than the chemical binding of 
the moderator. 

Graphite as a moderator is ideal for this purpose because the carbon atom is 
bound to the graphite lattice with an energy of about 20 ev. In all other respeats 
graphite oonforms well to the limitations imposed by the slowing down theory. 

A monoenergetic source of neutrons whose energy is between 20 ev and a lower 

limit of 1 ev imposed by restrictions to be explained later* was not obtainable . 

A monoenergetic source of neutrons of about 24 kev was used* therefore, the spatial 

distribution of 4.8 + 1 ev^^ gold resonance neutrons was obtained and compared by 

the Fermi theory to the 1.44 ev indium resonance neutron spatial distribution obtained 

(5) 

by Roberts, Hill and McCammon. 7 
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II Theory 
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It may be of interest to review the assumptions and limitations imposed in the 
derivation of the slowing down theory. Consider a fast neutron that makes an elastia 
•o llis ion with a stationary unbound nucleus and transfers a portion of its energy 
to the nualeus. The energy lost by the neutron depends upon its initial energy, the 
angle through whiah it is deflected and the mass of the nualeus. If the target atom 
is heavy, aarbon for example, and the kinetia energy assoaiated with it aonsists of 
t&'e thermal energy, its veloaity aan be negleated when aompared to the veloeity of 
a 1 ev neutron.^ Further, if in a center of mass system the probability of the 
neutron being scattered in any small solid angle is a constant, it may be shown 

that the average loss of energy by the neutron depends only on its initial energy and 
the mass of the target nucleus. The average loss of logarithm of the energy per col- 
lision, ^ » is independent of the initial neutron energy and is defined as amct'tgr, 
equal tos 



/ 
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where 




and 

s initial energy of the neutron 
E£ = energy of the neutron after collision 
A = mass of the target nucleus. 



( 1 ) 



When fast neutrons are released at a constant rate from a point in an infinite 
medium in which the above described collisions occur, the neutrons will diffuse away 
from the source and slow down until they are in thermal equilibrium with the medium. 

If the slewing down is treated a 3 a continuous process and the number of neutrons lost 
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by absorption during slowing down is negligible, the number of neutrons that slow 

•Z 

down past energy E in one second per cm at a distance r from the source is 
given by 



q(r,E) 



^ jo 

[Wr'( E J } 3 / 2 *’ ^ (3) 



( 2 ) 



Q 0 « source strength number of neutrons per seo. 

q(r,E) = slowing down density i.e, number of neutrons slowing down past 

energy E per sec cm^ 



T(E) s Fermi age 




(1 

s 





( 3 ) 



N 






Ei 



E. 



if cr s is 



« N(T^ 83 macroscopic cross section for scattering 

«* number of target nuclei per cm 
“ microscopic cross section for scattering 

■ initial neutron energy 
“ energy of neutron at detection. 

independent of neutron energy, 'C (E) may be integrated to give 



'C(B) 






(4) 






Treating the slowing down of neutrons as a continuous process rather than a dis- 
continuous process introduces the largest error when one attempts to find the slowing 

down density for age zero or at distances far from the source for any energy above 

, (2,7) 

thermal energies o 

Because the slowing down density is related to the neutron flux, foil activation 
measurements can be used to determine the slowing down density. The relation is 
given below* 



Q(r,E)dE - . 3^1 



(5) 



where Q(r,E) denotes the neutron flux and is identically equal 
to n(r,E)v 

n(r,E) * number of neutrons per cm® in energy range dE 
v ® scaler value of the neutron velocity cm per sec. 



For the purposes of this experiment, foil activations determined at various 
distances from the source were plotted, and (E) was determined from the slope. 

The 4.8 ¥ 1 ev neutron age can be compared to the age of 147 cra^ for 1.44 ev neutrons 
determined by Roberts, Hill and McCamraon by equation (4) 
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Ill Equipment and Techniques 

Except for the gold foils and counting instruments the equipment was the same 

as used by Roberts* Hill and McCammon and much of the following data were obtained from 

(5) 

their report 

lo Graphite column . 

The 6* x 6* x 11 • column shown in Fig. (I) was composed of 4” x 4" 
graphite stringers* the central stringers along the vertical axis were 
removable for inserting the foils* Its size was such that it was effective- 
ly infinite in size for the neutron ages investigated* The average density 
of the graphite was 1.652 g/cm^. The scattering cross section of the carbon 

• a co in " 24 2 

is 4.59 x 10 cm „ 

2* Neutron source. 

Practically monoenergetic neutrons were obtained from a source con- 
sisting of about 1 curie of antimony inside of a spherical shell of beryl- 
lium. The shell had an overall diameter of about 2.5 cm and weighed 9.24 
grams. 

Because the antimony has a rather high capture cross section for neutrons* 
a depression of the flux might be expected within a couple of mean free paths 
from the source. The absorptions of its own neutrons may also affect the 
half life of the antimony. This change is shown in thje appendix to be less 
than 0.0067^ and will therefore not affect the 60 day half life of the source 
appreciably. 
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Detecting Foils o 

Gold foils measuring 4 x 6.35 cm and 103.1 rng/crn^ thick were used to 
detect the 4.8 ev neutron flux. The maximum difference in the weight was 
6/S. During the exposure* the foils were enclosed in cadmium boxes shown in 
Fig.(ll) to shield them from thermal flux. 
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Two thin wall glass Geiger Mueller tubes enclosed in a 2" thick lead chamber 
to reduce the background counts were connected through preamplifiers to scales of 64. 
The tubes were normalized with activated gold foils over a range of counting rates 
that varied between 5,000 to 15,000 counts per minute. The normalization factor was 
found to be constant over this range of counting rates. 

Frequent counts were taken of two long half life foil standards to ascertain 
that the counting efficiencies of the instruments were constant. The two standards 
used were a Ru"* - ^ ft emitter with a half life of one year and a natural uranium 



foil. 
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Neutron Flux Measurements 

The isotopic neutron flux was detected in the following manner* Foils protected 
from the thermal flux by cadmium boxes were inserted at a number of positions along 
the principal axis of the graphite column* During the exposure time which consisted 
of about three days, the foils were oriented in the plane perpendicular to the 
principal axis of the column* After this aotivation time they were removed from 
the oolumn and from the cadmium boxes, wrapped inside of aluminum cylinders, shown 
in Fig* III and position on the Geiger Mueller tubes by a jig shown in Fig. IV, ? 
and counted for about eight hours* Immediately after one side of the foils had 
been counted, they were inverted and a similar count taken from the other side* 

By use of the normalization equation given in appendix 3 the counts obtained from 
the foils were converted to activity based on a 21*65 day exposure* The calculated 
activity of both sides of the foil was averaged to eliminate the current component 
of the flux* It is this averaged activity that is proportional to the isotropic 

( 5 

flux' * '* These data are given in column 7 of Table I* 

To evaluate the condition of the graphite column, indium foils were similarly 
exposed at a few positions in the column and counted* The calculated activities of 
the foils were compared with the activities reported by Roberts et al for these 
positions* A tabulation of these measurements are given with their data in Table 
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TABLE II 






Distance of foil from 



Measured saturated activity of indium 
foil covered with, cadmium 



neutron source 


Check Points 


\ 

Roberts, Hill, 


* 

> McC amnion 


r-cm 


2 2 

r -cm 


Act(r,E) x 1.0912 


Act(r,E) 


LogioAct(r,E) 


3.551 


12.61 


3590 


3529 


3.54720 


6.091 


37.10 


3236 


3162 


3.50001 


8.631 


74.49 


2924 


2836 


3.45271 


n.171 


124.79 




2563 


3.40870 


13.71 


187.96 




2240 


3.35017 


16.25 


264.06 




1906 


3.28007 


18.79 


353.06 




1581 


3.19846 


21.33 


454.97 




1276 


3.10584 


23.02 


529.9 




1122.8 


3.05030 


23.82 


567.4 


1030 


1051.2 


3.02168 


28.10 


789.6 




7H.3 


2.85205 


28.90 


835.20 




656.3 


2.81713 


33.18 


1100.9 




410.8 


2.61363 


33.98 


1154.6 


365.0 


373.8 


2.57262 


38.26 


1463.8 




217.O 


2.33646 


39.04 


1524.1 




197.4 


2.29542 


43.34 


1878.4 




113.3 


2.05423 


48.42 


2344.5 




52.09 


1.71675 


49.27 


2427.5 




47.67 


1.67822 


53.50 


2862.3 




26.65 


1.42570 


54.35 


2953.9 




21.50 


1.33245 
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Experimental Results 



A plot of the logarithm of the foil activity versus radii squared is 
given in Fig. V. A Gaussian distribution, equation 2, in these coordinates 
will appear as a straight line vith a slope of - 1/4 T( E ) . 

The neutron age may he obtained quite reasonably as 128 t 2 om c from the 
slope of the central portion of the curve. Near the origin much of the devi- 
ation from equation 2 may be attributed to the higher energy absorption 

5 

resonances in gold . At a radius of about 22 cm this effect becomes 
negligible and the curve conforms quite veil to a simple Gaussian out to a 
radius of about 45 cm. 

The experimental difference between the indium resonance and gold 
resonance neutron ages is 19 t 2 cm 2 . By the use of equation 6 the calculated 
age difference between 4.8 t .1 ev neutrons and 1.44 ev neutrons is 18.6 - .3 
cm 2 . These results then indicate that within the accuracy of this experiment 
the Fermi slowing down model is valid for neutrons slowing down in graphite 
from 4.8 ev to 1.44 ev. 



Corrections . 

The use of a spherical neutron source and finite size detector foils instead 

of a point source and point detectors have no effect on the results of this ex- 

(5) 

periment. Roberts , et al have shown that the effect of the spherical neutron 
source is to change the constant of proportionality between the foil activation 
and the neutron slowing down density. They also stated that a similar analysis* 
which is given in Appendix 2, would demonstrate that the only effect of the finite 
foils is to again change the proportionality constant. Because the absolute value 
of the neutron slowing down density is not of importance in this experiment, the 

' Q r ir‘ 1 1 V. i - ■ i J. *<~~J J.‘ ! 










corrected constant was not evaluatedo 

Unfortunately foils become activated by neutrons whose energies are not with- 
in the resonance energy range, i.e. resonance energy minus to resonance energy plus 
one half the average energy loss per collision* To correct for this non resonance 
activation it is required that the variation of capture cross section with neutron 
energy be known. 

Such a correction was attempted following the analysis given in Appendix 4 
and using for the capture cross section for gold the Breit Wigner formula 





cr 0 p 2 

£(e - E a ) Z - pz/ij 



where E 0 ** 

r * 

C 71 - 



4.8 ev, the neutron energy for resonance peak 

.481 ev, the nuclear energy level width 

2700 barns, the capture cross section at resonance peak. 



This procedure resulted in a corrected age for the gold resonance 

neutrons of 127 cm*”. Because the values of /~7 is perhaps much smaller than 

2 2 

o481 ey, the corrected age should perhaps be somewhere between 127 cm and 128 cm • 
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APPENDIX 



1* To show that the absorption of neutrons affects the antimony half life a 
negligible amount, the assumption is made that every neutron produced by the 
source is self absorbed* Wattenberg^ states that approximately one half of 
the antimony disintegrations produce the 1*8 Mev y ray* In the same article the 
cross section of the y (Be-n) reaction is given as 9*7 x 10“^® cm^* 

The approximate rate of forming neutrons is given by 



nA /, ~ N Be a ”aBe T Be,. 
~~ 2 — '• L " 9 ) 

Tg e “ thickness of Be shell 
A 8 decay constant for antimony 
N A s | rate of decay of the antimony | 
s *317 x 10“^ Be atoms per cm? 



N 



Be 



* 0*4365 cm, thickness of Be 



rate - M A 
2 



-0*000134. 

" ® ) s 0.000067 N A. 



The actual rate of decay of the antimony may be obtained from 



a - An + (1 - e “0*°00134^ 

-f|- = - A N(1 - 0.000067) 

and 

A apparent * (1 - 0*000067) A 

2* To show that no foil geometry corrections are necessary, the following assumptions 
are made* 



a* Point neutron source, proof that the actual source appears as a point source 
is given in 0RNL-201* 



b* Circular, thin foils« 

c« The spatial resonance flux distribution is Gaussian. 
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Since and 'V (E 0 ) are to be constant the total foil activation varies with. 

-r 2 /4U (E ) 

distance from the source as e 0 • 

3* All foils were corrected to the activation expected if they had been placed in 
the column the first day, exposed for 21o9 days and counted immediately* 

N a « number per cm 3 of Au^ -97 
R * number per cm^ of active gold atoms 

A^ s antimony decay constant m 0*0155/day* 

A g ■ gold decay constant “ 0*2567/day* 

The source strength varies with time as 

* Ayfc 

If (t) = A e 1 

TQO T_Q 7 

If the number of Au atoms formed is small compared to the number of Au A9 ' 
originally present,the rate of formation of Au^ 99 may be expressed as 

= - A 2 N+-N a O- a f(t) 



dH 

dt 



but N a <r^ is considered constant rfcheref ore let A N a o~^ s B 



The boundary condition is that N s 0 when t * 0* 

The number of Au nuclei formed at any time t, is then 

-Aiti -Ag^ 



B 



N i s rvy 

A 2 Ai 



(® 



— a 



) 



When the foils are removed from the column »the gold decays as 

dN \ 

— = -H A, 



r- rfro 

m — ■-» f T“ 



■Xl; iL k.'J Jl X'OiV — 









KJ x J. 






The number present Ng at the end of this cooling interval tg will be 



N, 



•^l t l "^2*1. “^2*2 

D (e - © ) © 



^2 " \ 



The counts obtained, C, in the counting interval tg will depend upon the counter 
efficiency d * 

t- t* 



4tS - 



- Agt 
e dt 



C = Of 



B 



(• 



- Aiti -A 2 ti m ^2^2 " < ^2' b 3 



- e 



) e (1 - e ) 



^2 “ ^1 

The maximum number of counts that could be obtained would occur if t^ were equal to 

jgn * z /\ 

“ V" 1 \ * *2 ®<l ua l bo zero and tg equal to infinity* For the decay constants 

A i - / 2 

given for antimony and gold, t^ for maximum activation in the column is 21*9 days* 



at n -21*9/}, -21*9/^ „ 

Vx ■ - ~ — r - (• 1 -> 2 ) 



A 2 -Aj 



By forming the ratio of maximum counts to the actual counts and multiplying by the 
gold decay constant the activity normalization factor for 21*9 days exposure time 
and zero cooling time is obtained* 



^ 2 ^ 



-21.9 A, -21.9 A, 



A * *1° * 111 *1 

max " ^2 t l " ^2 t 2 “A 2*3 

(e - e ) e (1 - e ) 






To correct the activity to that which would have been attained if the foils had 



been inserted into the column on the 1st day, Act mRY must be multiplied by the 
" Aj^t^ 

source factor l/e ® The is the time elapsed between the insertion of the 

first foil on the first day and the insertion into the column of the foil one 
wishes to normalize® 

The normalizing equation becomes 

2l®9 -21®9 



Act 



>5 0 



(max 






- e 



j 



“Ait 



1st day) (©” *ri m e " a 2“1) 
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4® To account for low energy gold activation the following analysis is applied® 

d-/2. s differential solid angle 

u “ cos 0 

ds s differential area 
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Assumptions j 

a* Isotropic flux* 

b* The scattering cross sections of the shield and foil materials can be 
neglected* 

c* Absorption of neutrons of energy E does not perturb the neutron flux at 
any lower energy* 
d* Plane geometry* 

Consider the foil at position r Q 

(1) y (E,r 0 )dE - number of neutrons with energies between E and 
E j- dE which cross unit area per second* 

(2) If (E,r Q )dE u a number of neutrons with energies between E 

and E +• dE which cross unit area normal to the foil per second 

in the solid angle d -A- s 2 ff du* 

-N_ O" z/u 

(3) If (E,r Q )dE u © s the number of these neutrons which 

penetrate the shield to a depth of z cm per second if dr^ is the 
absorption cross section and N^ is the number of absorbing atoms / cm® 
of shielding material* 

(4) Of the neutrons that penetrate this far the number absorbed in 
layer dz cm thick per second is 

-NjCTjz/u 
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or as 2 if du s d^ 
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du 
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(5) f(E, r 0 )dBu 



du 



z /u 



( -J!f dz) 



r -JSL<r T /u 7 

s (E,r 0 ) u dE - ® J s 



number per see 

absorbed in thickness in angle du* 

(6) ^(E,r 0 )u — dE e ^ ^ ^ jj number of neutrons with energies 

between E and E -f- dE getting through one cm^ of shield per second 
with angle du about u* 

(7) By applying similar reasoning to the neutron absorption in the foil 
and multiplying the result by a factor of 2 to account for aotivation 
from the other side the activation becomes 



-N CT l /u f -H *cr *T */u 7 

y(E»r 0 )u © 1 1 1 p - © J 



dE du 



where N'o-’T* refers to the gold foil* 

(8) Since (j> (E,r 0 ) is independent of u, the equation in (7) may b© 
integrated over u with limits of 0 to 1* By a change of variable 
u s -y- the activation becomes 
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■where 
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dt 



CT^ and cr* are functions of energy 
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^ and £ to be evaluated for graphite* 
s 

(9) The correcting factor used is 

-,(4.8 0.34) ev 
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